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1 Introduction 

Deeply virtual Compton scattering (DVCS), 'j*{qi)p{Pi) — > 7(q'2) ^(-Pj), is viewed 
as the cleanest process to access generalized parton distributions (GPDs), which 
encode a partonic description of the nucleon, cf . Refs. [21 H] . In the kinematics of 
HI and ZEUS collider experiments at HERA, the DVCS cross section is to a large 
extent dominated by the flavor singlet part of the helicity conserved Compton form 
factor (CFF) ^: 



'iiW,t,Q^)^'^e\^iU = A\Q^ 



(V 



•^Talks given by K.K. at DIS 2008, 7-11 April 2008, London and D.M. at International 
Workshop on Diffraction in High-Energy Physics, 9-14 September 2008, La Londe-les-Maures, 



France [5]. 
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Here, W is the c.o.m. energy, A = P2 — -Pi is the momentum transfer, — = 
is the incident photon virtuahty, and ^ ~ is a Bjorken-hke scahng variable. 

The CFF ^ factorizes further into a convolution of the partonic, i.e., hard 
scattering, amplitude C = (^C, *^C) and GPDs H = (^if, ^if), (E=singlet quark, 
G=gluon), 

^(e, t, Q') = fdx C(x, QVf^', H{x, v = ^,t, i?) , (2) 

where the skewness parameter rj = —A ■ q/{Pi + P2) ■ Q is set equal to ^. The 
factorization scale fi separates short- and long-distance dynamics and is often 
taken as /i = Q. The scale dependence is governed by evolution equations. They 
also tell us that the evolution effect depends on the (x, rj) GPD shape itself. Note 
that gluons do not directly enter the DVCS amplitude at leading order (LO), but 
rather drive the evolution of singlet quarks. 

Since the momentum fraction x is integrated out in the amplitude ([2]), GPDs 
cannot be directly revealed. In the quest for a realistic GPD model, GPD values 
along the two trajectories rj = and rj = x play a prominent ro 1(0. Namely, GPDs 
on these trajectories are given at LO by DIS structure function and imaginary part 
of DVCS amplitude, respectively. In this sense they are experimentally measurable. 
Thus, a modelling strategy that places emphasis on these trajectories has a good 
chance to be efficient in a global DVCS fit and to capture the physical GPD 
content, giving us a link to a partonic interpretation. 

It is the objective of this work to find GPDs that satisfy the well-known the- 
oretical constraints [3l H] and provide a good fit to all available HI and ZEUS 
DVCS data. Our particular concern is the LO description, where so far this goal 
has not been reached. Since the real part of the CFF 7i can be calculated from a 
dispersion relation, our outcome might be relevant for the real part in fixed target 
kinematics, too. Beyond LO these DVCS data can be described within specific 
GPD models, see, e.g., Refs. [6l[7j. 



2 Conformal moment- space approach 

It is convenient to work with conformal GPD moments. For integral conformal 
spin j + 2 they are defined by convolution with Gegenbauer polynomials Cj'(x), 
e.g., for quarks: 

H]{r^, t, f,') ^ ^M^Wl Cf{x/v) H^ix, V, t, /i^) , (3) 

■' 2-'+lr(j+3/2) J-1 ■' 

^ These trajectories appear in a GPD sum rule family, which can serve as powerful modeUing 
tool 0. 
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where q is the flavor index. The normahzation ensures that in the forward hmit 
A — > the conformal moments simply reduce to famihar Melhn moments of parton 
distribution functions (PDFs). Some of the advantages of working with conformal 
moments are 

• Conformal moments evolve autonomously at LO and, in a special factoriza- 
tion scheme, even at next-to-leading order (NLO). 

• Powerful analytic methods of complex j plane are available (similar to com- 
plex angular momentum). This also allows to build a stable and fast routine 
for fitting. 

• New possibilities for GPD modelling, which makes direct contact to the t- 
channel SO (3) partial wave expansion, Regge phenomenology, and lattice 
measurements. 

In moment space the convolution formula ^ yields formally a divergent series 
over the integral conformal spin. Analogously to a SO (3) partial wave expansion, 
it can be resummmed by means of a Mellin-Barnes integral representation [H [9|, [7] 



1 pc+ioo 
Z/l J c—ioo 



i + tan ( ^ 



C,(QV/i2,a,(/i))if,(e,t,/i') 



(4) 

GPD moments Hj, cf. Eq. ([3]), can be decomposed in t-channel SO (3) partial 
waves, i.e., Wigner matrices (iQ^(cos^), which are labelled by angular momentum 
J and hadron helicity differences {y = 0, ±1) [IHIEIE]- The cosine of the scattering 
angle 9 might be approximated by —l/r]. We include an effective leading Regge 
pole at a{t) = ao + a't in the partial wave amplitudes and parameterize their 
residual t-dependence by a dipole with cut-off mass Mj (or, alternatively, with an 
exponential t-dependence). Our model (for assignments see Ref. p]) is given at a 
input scale Qq and reads for integral j: 

j=j™ ^y^) (l 



where hj are strengths of partial waves. The PDF Mellin moments hj ''^^ /{j + 
1 — cto) are fixed from DIS, the remaining parameters are constrained by DVCS 
fits. 

In practice the parameter space in the ansatz ([5]) must be further reduced. We 
employ three GPD models: i. taking only the leading J = j ' + 1 SO (3) partial wave 
(1-PW), ii. including the next-leading J = j — 1 one (nl-PW), and iii. performing a 
model dependent resummation of SO (3) partial waves (S-PW). To set up the third 
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Figure 1: Fit of S-PW model at LO to DVCS and DIS data from [El [H [ig [16] . 



model, we assume that the GPDs in the vicinity of x = 77 behave as ^(f^)^""'-*-'- 
Relevant for us is the small r] behavior of the resulting conformal moments for 
complex- valued j which is given by a series of 'conformal' daughter poles 
j = a{t) — 1 — n with n = 2, 4, 6, ■ ■ ■ . All of them contribute to the leading Regge 
behavior of the CFF Ti. The strength of the non-leading SO (3) partial waves 
in the nl- and S-PW model for sea quarks and gluons is controlled by skewness 
parameters. These parameters allow us to adjust the normalization of CFFs. Note 
that Mellin-Barnes and "dual" [lOJ GPD parameterizations are related. In the 
latter p = j + 1 — J is taken as a parameter and a 'forward-like' momentum fraction 
^-integral replaces the j-integral (jl]). Our 1- and nl-PW models are equivalent to 
a minimalist (p = 0) and a minimal (p = {0,2}) "dual" model, respectively [TT]. 



3 Results 

The 1-PW model is the simplest one with a rigid normalization, i.e., Ti. at Qo s'^nd 
t = is given by the unpolarized sea quark PDF enhanced by the Clebsch-Gordan 
coefficient 

. 1.5. (6) 
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Figure 2: (a) the skewness ratio r, cf. Eq. (jTj), for fixed x = 10^'^ (thin) and 
alternative ratio R, cf. Eq. (IHl), for fixed W = S2 GeV (thick), (b) and (c) effective 
exponential t-slope (E]) versus Q and W , respectively, together with experimental 
data from Ref. [1^; (d) transverse width ( fTOl) of sea quarks at the input scale 
= 4GeV^ (thin) and evolved at = 10 GeV^ (thick). The model parameter 
for 1 (dotted), nl (dashed), and S (sohd) S0(3)-PW models, see text, are obtained 
from LO fits such as in Fig. [H 



Within such a simplified model the t-slope is poorly described at LO, since it is 
used to adjust the normalization, whereas at NLO or beyond good fits are obtained 
[7]. To have a LO description, one has to include non-leading S0(3) partial waves. 
Both of our flexible GPD models (nl-PW and S-PW) allow for good DVCS fits at 
LO, illustrated for the S-PW model in Fig. [1], and beyond. The t-dependence is 
well described with a' = 0.15/GeV^ and Mj ~ 0.6 GeV, taken at the initial scale 
Qo = 2 GeV. Note that to some extent these parameters are correlated. The last 
panel shows fit to DIS structure function F2, which fixes the normalization /i^'*'^ 
of leading S0(3) partial wave and intercept ao- We add that gluonic GPD is not 
well constrained, e.g., its t-slope can be taken from J/ip data ^12j . 

We define a skewness ratio of singlet quark GPD on the r] = x and r] = 
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trajectories: 

Resulting from our LO fit in Fig. [1], r versus is plotted in Fig. [2](a) for the 
nl- and S-PW model as thin dashed and solid curves, respectively. One notes 
that for both models r ~ 1 over the wide lever arm, i.e., a skewness effect is 
almost absent at LO and the models are not distinguishable. In these models the 
gluonic r ratio is considerably smaller than one and can even be negative. The 
ratio R (thick curves) of imaginary parts of DVCS amplitude and DIS structure 
function might be considered as an observable [13]. This alternative measure of 
the skewness effect in the LO approximation reads 

^ H{x,x,t = 0,Q^) 



H{2x,0,t = 0,Q^)' ' ' 

Note that R ^ 2"(°'2')r with an intercept a(0, ^ 4 GeV^) ~ 1.2. Numerous 
GPD models and the claim r ~ 1.5 [T7], which is implemented in the 1-PW model 
(dotted curves) and arises from the Clebsch-Gordan coefficient IQ, are disfavored 
by our LO fits. As mentioned above, we are even concerned about their uses in 
fixed target kinematics, if described at LO [18j. We emphasize that 1-PW models 
with r ~ 1.5 can describe the DVCS data beyond LO [7], where the corresponding 
gluonic ratio is about one. 

The t-dependence of the DVCS data can be directly fitted within an exponential 
ansatz, where a Q dependent t-slope and a vanishing a' has been found We 
display in Fig.[2](b) and (c) the HI measurements [13] and the effective exponential 
t-slope 

^eff ^ 1 1^ ^(Iy,t = -0.8GeV^62) 
-0.7 GeV^ ^^^Ep(iy,t = -O.lGeV^ Q2) ' 

evaluated from a E-PW model (solid) with dipole ansatz and a' = 0.15/GeV^ 
as well as a nl-PW model (dashed) with exponential ansatz and a' = 0. Models 
describe the present DVCS data set with x^/^-O-f- ~ 101. /98 and x^/d.o.f. ~ 
98. /98, respectively. We emphasize that the differences between the nl- and E-PW 
model, visible in Fig. [21 mainly originate from the implementation of t-dependence. 
The decrease of the t-slope with growing photon virtuality Q, shown in panel (b), 
entirely arises from evolution. The flatness of the t-slope with respect to the W- 
dependence, see panel (c), does not exclude a small a' value. Note that due to the 
evolution the a' value at the scale = 10 GeV^ is smaller than the input value 
a' = 0.15/GeV^. We add that from our fits a large value of a' ~ 0.8/GeV^ for the 
dominant sea quark contribution is disfavored, cf. Ref. 



6 



Present small-xej DVCS data are well described by different GPD models and 
thus a partonic interpretation w.r.t. the skewness dependence, i.e., the distribution 
of SO (3) partial waves, and the t-dependence can not be unique. For the transverse 
width [I9j 

(6^) (x, Q2) = 4^ In H{x, r/ = 0, t, (iq) 

of sea quarks this is illustrated in Fig. [2] (d). Here the differences for the transverse 
width of the dipole (solid) and exponential (dashed) t-dependent model are caused 
by their extrapolation descriptions of experimental data to t = 0. We also show 
the evolution of the transverse width from = 4GeV^ (thin) to = lOGeV^ 
(thick). Thereby, the slope a' decreases, compare thick and thin solid curves. 

We remark that our models also reproduce the preliminary HI result on the 
dominant twist-two cos(j) harmonic for the beam charge asymmetry [20j. Since this 
asymmetry is proportional to the real part of a OFF combination, its description 
gives us confidence in the Regge behavior of the DVCS amplitude. Unfortunately, 
we could not obtain a useful bound for the OFF S or the related anomalous 
gravitomagnetic moment, which would shed light on the partonic decomposition 
of the proton spin. 

In conclusion, we successfully describe the small xbj DVCS data within flexible 
GPD models from the perspective of t-channel physics or Regge phenomenology. 
Thereby, a LO description requires the inclusion of non-leading SO (3) partial 
waves, while beyond this order they are not essential. The interplay of both t 
and ^ dependence with evolution is in our GPD models compatible with data. 
In a partonic LO interpretation the dominant sea quark GPD possesses almost 
no skewness effect over a wide lever arm. This feature can be realized within 
various models. We also found that gluons, driving the evolution at small x, are 
relatively suppressed (even a negative gluonic GPD on the trajectory t] = x). The 
functional form of the t-dependence in the GPD models can not be fully pinned 
down from present DVCS data. In particular, a small, however, non zero, a' can 
not be excluded and the residual t-dependence might be given by a dipole ansatz. 
For a global fit of DVCS data, a good LO description of small-xsj DVCS data is 
valuable at present. Whether one prefers a quark interpretation of DVCS data in 
a specific factorization scheme, in which absence of skewness effect holds at any 
order, or likes to resolve the gluonic content within radiative corrections in the 
standard scheme, is primarily a convention, yielding a GPD reparameterization. 
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